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Abstract 

A measurement of angular correlations in Drell-Yan lepton pairs via the 0* observable is presented. 
This variable probes the same physics as the Z/j* boson transverse momentum with a better experi- 
mental resolution. The Z/7* ->■ e+e~ and Z/7* ->• decays produced in proton-proton collisions 
at a centre-of-mass energy of y/s = 7 TeV are used. The data were collected with the ATLAS detec- 
tor at the LHC and correspond to an integrated luminosity of 4.6 fb _1 . Normalised differential cross 
sections as a function of <f>* are measured separately for electron and muon decay channels. These 
channels are then combined for improved accuracy. The cross section is also measured double dif- 
ferentially as a function of <f>* for three independent bins of the Z boson rapidity. The results are 
compared to QCD calculations and to predictions from different Monte Carlo event generators. The 
data are reasonably well described, in all measured Z boson rapidity regions, by resummed QCD 
predictions combined with fixed-order perturbative QCD calculations. Some of the Monte Carlo event 
generators are also able to describe the data. The measurement precision is typically better by one 
order of magnitude than present theoretical uncertainties. 
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Abstract 

A measurement of angular correlations in Drell-Yan lepton pairs via the 4>* observable is presented. This variable probes 
the same physics as the Zj^j* boson transverse momentum with a better experimental resolution. The Z/^y* e + e~ 
and Z/j* — > decays produced in proton-proton collisions at a centre-of-mass energy of y/s = 7 TeV are used. 

The data were collected with the ATLAS detector at the LHC and correspond to an integrated luminosity of 4.6 fb . 
Normalised differential cross sections as a function of (j>* are measured separately for electron and muon decay channels. 
These channels are then combined for improved accuracy. The cross section is also measured double differentially as a 
function of </>* for three independent bins of the Z boson rapidity. The results are compared to QCD calculations and 
to predictions from different Monte Carlo event generators. The data are reasonably well described, in all measured Z 
boson rapidity regions, by rcsummcd QCD predictions combined with fixed-order perturbativc QCD calculations. Some 
of the Monte Carlo event generators are also able to describe the data. The measurement precision is typically better 
by one order of magnitude than present theoretical uncertainties. 
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1. Introduction 

In hadron collisions at TeV energies the vector bosons 
W and Z/j* are copiously produced with non-zero mo- 
mentum transverse to the beam direction (pt) because 
of radiation of quarks and gluons from the initial-state 
partons. In this context the signatures Z/j* -4 e + e~ 
and Z/j* — > provide an ideal testing ground for 

QCD due to the absence of colour flow between the initial 
and final state [1-3]. The study of the low p^ spectrum 
(pf < mz), which dominates the cross section, has im- 
portant implications on the understanding of Higgs boson 
production since the transverse-momentum resummation 
formalism required to describe the Z/j* boson cross sec- 
tion is valid also for the Higgs boson [4-7]. A precise un- 
derstanding of the spectrum is also necessary to further 
improve the modelling of W boson production in QCD cal- 
culations and Monte Carlo (MC) event generators, since 
the measurement of the W mass is directly affected by 
uncertainties in the pt^ shape [8, 9]. 

The transverse momentum spectra of W and Z/j* 
bosons produced via the Drell-Yan mechanism have been 
extensively studied by the Tevatron collaborations [10-14] 
and, recently, also by the LHC experiments [15-17]. How- 
ever, the precision of direct measurements of the Z/j* 
spectrum at low p^ at the LHC and the Tevatron is lim- 
ited by the experimental resolution and systematic uncer- 
tainties rather than by the size of the available data sam- 
ples. This limitation affects the choice of bin widths and 
the ultimate precision of the p^ spectrum. In recent years, 



additional observables with better experimental resolution 
and smaller sensitivity to experimental systematic uncer- 
tainties have been investigated [18-21]. The optimal ex- 
perimental observable to probe the low-p^ domain of Z/j* 
production was found to be 4>* which is defined [20] as: 

4>* n = tan(> ac0 p/2) • sin(0*) , (1) 

where acop = ti — A<fi, A(f> being the azimuthal open- 
ing angle between the two leptons, and the angle 9* is 
a measure of the scattering angle of the leptons with re- 
spect to the proton beam direction in the rest frame of 
the dilcpton system. The angle 9* is defined [20] by 
cos(#*) = tanh[(?/ _ — rj + )/2] where rf~ and ?y + are the 
pseudorapidities 1 of the negatively and positively charged 
lepton, respectively. Therefore, <fi* depends exclusively on 
the directions of the two lepton tracks, which are better 
measured than their momenta. The <j>* variable is posi- 
tive by definition. It is correlated to the quantity p^/mu, 
where mu is the invariant mass of the lepton pair, and 
therefore probes the same physics as the transverse mo- 
mentum p% [22]. Values of <fi* ranging from to 1 probe 



1 ATLAS uses a right-handed coordinate system with its origin at 
the nominal pp interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the 
IP to the centre of the LHC ring, and the j/-axis points upward. 
Cylindrical coordinates (r, <j>) are used in the transverse plane, </> 
being the azimuthal angle around the beam pipe. The pseudorapidity 
is defined in terms of the polar angle 8 as n = — lntan(6>/2) and the 
rapidity is defined as y = ln[(_E + p z )/(E —p z )]/2. 
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the distribution mainly up to ~ 100 GcV. The (f>* distri- 
bution of Z/-f* bosons has been measured in three bins of 
the Z boson rapidity (yz) by the D0 Collaboration using 
7.3 ffcr 1 of pp collisions at y/s = 1.96 TeV [23]. 

This Letter presents a measurement of the normalised 
tf>* distribution in bins of the Z boson rapidity yz using 
4.6 fb _1 oipp interactions collected at i/s = 7 TeV in 2011 
by the ATLAS detector. The normalised differential cross 
section is measured in both the electron and muon chan- 
nels in the fiducial lepton acceptance defined by the lepton 
(£ = e, /i) transverse momentum p^ > 20 GeV, the lepton 
pseudorapidity |t/| < 2.4 and the invariant mass of the 
lepton pair 66 GeV < m« < 116 GeV. Correction factors 
allowing the extrapolation of the cross section from the 
fiducial letpton acceptance to the full lepton acceptance, 
restricted to 66 GeV < m« < 116 GeV, are also presented. 
The reconstructed </>* distribution, after background sub- 
traction, is corrected for all detector effects. The measure- 
ments are reported with respect to three distinct reference 
points at particle level regarding QED final-state radiation 
(FSR) corrections. The true dilepton mass mu and </>* are 
defined by the final-state leptons after QED FSR ("bare" 
leptons), or by recombining them with radiated photons 
within a cone of AR = ^(Ar/) 2 + (A(f>) 2 =0.1 ("dressed" 
leptons), or by the final-state leptons before QED FSR 
( "Born leptons" ) . The Born definition corresponds to the 
full correction for QED FSR effects, so that it can be used 
for the combination of the electron and muon channels. 
The combination of the electron and muon channels is 
compared to QCD predictions obtained by matching re- 
summed and fixed order QCD calculations, as well as to 
the predictions of MC event generators implementing a 
parton shower (PS) algorithm. 

2. QCD predictions 

Non-zero p^* is mainly generated through the emission 
of partons in the initial state. In the high p^ region 
(Pt ^ m z) the spectrum is determined primarily by hard 
parton emission. Perturbative QCD calculations, based 
on the truncation of the perturbative series at a fixed or- 
der in a s , are theoretically justified and provide reliable 
predictions. The inclusive cross-section prediction is finite 
but the differential cross section diverges as approaches 
zero. In this limit (p^ <C mz) the convergence of the 
fixed-order expansion is spoiled by the presence of powers 
of large logarithmic terms which have to be resummed to 
restore the convergence. 

Differential cross sections calculated to 0(a 2 s ) are avail- 
able for Z/j* production through the Fewz [24, 25] and 
Dynnlo [26, 27] programs. The ResBos [28-30] gener- 
ator resums the leading contributions up to next-to-next- 
to-leading logarithms (NNLL) and matches the result to 
fixed-order calculations at 0(a s ). This is corrected to 
0(a 2 s ) using a k- factor depending onpf and yz [31]. In ad- 
dition, the ResBos generator includes a non-perturbative 
form factor that needs to be determined from data [32]. 



A slightly different approach has been proposed recently 
to describe the Tevatron Run II data by matching NNLL 
accuracy to MCFM calculations [33], with no apparent 
need for non-perturbative contributions [22, 34]. 

Similarly to resummed calculations, PS algorithms such 
as those used in Pythia [35] and Herwig [36] provide an 
all-order approximation of parton radiation in the soft and 
collinear region through the iterative splitting and radia- 
tion of partons. The Powheg [37-40] and Mc@nlo [41] 
event generators combine next-to-leading order (NLO) 
QCD matrix elements with a PS algorithm to produce 
differential cross-section predictions that are finite for all 
pfL The Alpgen [42] and Sherpa [43] event generators 
implement tree-level matrix elements for the generation of 
multiple hard partons in association with the weak boson. 
They are matched to parton showers either by a PS algo- 
rithm using re-weighting procedures [44, 45] or through a 
veto [42], in order to avoid the double counting of QCD 
emissions in the matrix element and the parton shower. 

3. The ATLAS detector 

The ATLAS detector [46] is a multi-purpose particle 
physics detector operating at one of the beam interaction 
points of the LHC. It covers nearly the entire solid an- 
gle around the collision region and consists of an inner 
tracking detector (inner detector or ID) surrounded by a 
thin superconducting solenoid providing a 2 T axial mag- 
netic field, electromagnetic and hadronic calorimeters, and 
a muon spectrometer (MS). 

Measurements in the ID arc performed with silicon pixel 
and microstrip detectors covering |?7| < 2.5. A straw- 
tube tracking detector follows radially and covers the 
range \r]\ < 2.0. The lead/liquid-argon electromagnetic 
calorimeter is divided into barrel (\rj\ < 1.5) and end- 
cap (1.4 < I77I < 3.2) sections. The hadronic calorime- 
ter is based on steel/scintillating tiles in the central re- 
gion (|7y| < 1.7), and is extended to |?y| = 4.9 by end- 
cap and forward calorimeters which use liquid argon. The 
MS comprises separate trigger and high-precision tracking 
chambers to measure the deflection of muons in a magnetic 
field generated by three large superconducting toroids ar- 
ranged with an eightfold azimuthal coil symmetry around 
the calorimeters. The high-precision chambers cover a 
range of \r/\ < 2.7. The muon trigger system covers the 
range |?7| < 2.4 with resistive plate chambers in the barrel, 
and thin gap chambers in the endcap regions. 

4. Event simulation 

MC simulations are used to calculate efficiencies and 
acceptances for the Z/j* — > £ + £~ signal processes and to 
unfold the measured spectrum for detector effects and 
for different levels of QED FSR. The Powheg MC gener- 
ator is used with CT10 [47] parton distribution functions 
(PDFs) to generate both the Z/7* e+e" and Z/7* 
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signal events. It is interfaced to Pythia 6.4 with 
the AUET2B-CTEQ6L1 tune [48] to simulate the par- 
ton shower and the underlying event. Generated events 
are re- weighted as a function of p^ to the predictions from 
ResBos, which describes the p^ spectrum more accu- 
rately [15]. Simulated events are also used to estimate 
background contributions. The electroweak background 
processes W — > iv and Z/-f* — > t + t~~ are generated using 
Pythia 6.4. The production of ti events is modelled using 
Mc@nlo and diboson processes are simulated using Her- 
wig. The event generators are interfaced to Photos [49] 
to simulate QED FSR for all of the simulated samples, ex- 
cept Sherpa which is interfaced to an implementation of 
the YFS algorithm [50, 51]. 

Multiple interactions per bunch crossing (pile-up) are 
accounted for by overlaying simulated minimum bias 
events. To match the observed instantaneous luminosity 
profile, the simulated events are re-weighted to yield the 
same distribution of the number of interactions per bunch 
crossing as measured in the data. The response of the 
ATLAS detector to the generated particles is modelled us- 
ing Geant4 [52], and the fully simulated events [53] are 
passed through the same reconstruction chain as the data. 
Simulated event samples are corrected for differences with 
respect to the data in the trigger efficiencies, lepton recon- 
struction and identification efficiencies as well as in energy 
(momentum) scale and resolution. The efficiencies are de- 
termined by using a tag-and-probe method similar to the 
one described in Section 4.3 of Ref. [54] based on recon- 
structed Z and W events, while the energy resolution and 
scale corrections are obtained from a fit to the observed Z 
boson line shape. 

5. Event reconstruction, selection and background 
estimation 

Events recorded during periods with stable beam condi- 
tions and passing detector and data-quality requirements 
are selected. At least one primary vertex reconstructed 
from at least three tracks is required in each event. 

Events in the electron channel are selected online by 
requiring a single electron candidate with a threshold in 
transverse momentum p^ that was increased during the 
data-taking from 20 GeV to 22 GeV in response to in- 
creased LHC luminosity. Electrons are reconstructed from 
a cluster of cells with significant energy deposits in the 
electromagnetic calorimeter matched to an inner detector 
track. Electron reconstruction uses track refitting with a 
Gaussian-sum filter to be less sensitive to bremsstrahlung 
losses and improve the estimates of the electron track pa- 
rameters [55, 56]. The typical angular resolutions in the 
electron direction measurements are 0.6 mrad for (j> and 
0.0012 for ?/. The highest and second highest pt electrons 
arc required to have a transverse momentum p^ > 25 GeV 
and pf> > 20 GeV, respectively. The electron pseudora- 
pidity must satisfy |?7 e | < 2.4 with the calorimeter bar- 
rel/endcap transition region 1.37 < \if\ < 1.52 excluded. 



Electrons arc required to pass "medium" identification cri- 
teria based on shower shape and track-quality variables, as 
described in Refs. [57, 58]. The criteria are re-optimised for 
both higher pile-up conditions and higher instantaneous 
luminosity in 2011. 

Events in the muon channel are selected online by a trig- 
ger requiring a single muon candidate with pt^ > 18 GeV. 
Muons are identified as tracks reconstructed in the muon 
spectrometer matched to tracks reconstructed in the in- 
ner detector and are required to have pi^ > 20 GeV and 
1 77^ | < 2.4. Only isolated muons are selected by re- 
quiring the scalar sum of the pt of the tracks within a 
cone Ai? = 0.2 around the muon to be less than 10% of 
the muon px- Muons are required to have a longitudi- 
nal impact parameter with respect to the primary vertex 
less than 10 mm to reduce contributions from cosmic-ray 
muons and in-time pile-up. In addition, the transverse im- 
pact parameter of the track with respect to the primary 
vertex divided by its uncertainty must be smaller than ten 
to reduce non-prompt muon backgrounds. The typical an- 
gular resolutions in the muon direction measurements are 
0.4 mrad for </> and 0.001 for rj. 

Z/7* — > £ + £~ events are selected by requiring two oppo- 
sitely charged same-flavour leptons with an invariant mass 
66 GeV < mu < 116 GeV. After these selection require- 
ments 1.22 • 10 6 dielectron and 1.69- 10 6 dimuon candidate 
events are found in data. 

Background contributions from Z/j* — > t + t~, W — > 
tv, ti and diboson production are estimated using MC 
simulations. The cross sections are normalised to next-to- 
next-to-leading-order (NNLO) predictions for Z/7* and 
W production using Fewz, NLL-NLO predictions for ti 
production [54] and NLO predictions for diboson produc- 
tion [59]. For both the e + e~ and channels, the main 
background at high <j)* values arises from ti and diboson 
production. 

At low 0* values the background is dominated by multi- 
jet production, where a jet is falsely identified as a primary 
e or [i. In this case the background is determined by data- 
driven methods. A data event sample dominated by jets 
faking electrons or muons in the final state is employed 
to determine the shape of the multi-jet background. For 
the e + e~ channel, the multi-jet sample is obtained from 
electrons failing the medium identification criteria. In or- 
der to assess systematic uncertainties in the shape of the 
multi-jet background, an alternative multi-jet control sam- 
ple was also selected using non-isolated electrons. For the 
H + pT channel, the multi-jet sample is extracted by in- 
verting the isolation requirement on muons. The uncer- 
tainty in its shape was studied by comparing same-sign 
and opposite-sign dimuon events. The normalisation of 
this multi-jet background template is determined by ad- 
justing the sum of it and other background and signal MC 
predictions to data as a function of the invariant mass 
spectrum of the dilcpton pair. An extended dilcpton mass 
range, 50 GeV < m u < 150 GeV (200 GeV for electrons), 
was employed to better constrain the off-resonance region 
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Table 1: The measured normalised differential cross section l/<r fid ■ d<r fid /d0* in bins of tf>* for Z/"/* — > e + e — and Z/y* — y A* + A t_ channels. 
The cross sections, which are to be multiplied for convenience by a factor /, are reported with respect to the three different treatments of QED 
final-state radiation. The relative statistical (5 s tat) an d total systematic (<5 S y S ) uncertainties are given in percent. The overall point-to-point 
uncorrelatcd additional uncertainty in QED FSR of 0.3% is not included. 







Z/ 7 * - 


-+ e+e 












Z/7* - 












< 


K 


l/(T fld 








<5stat 


6 


sys 


l/cr fid 


d ( r fld /d^; 






<5stat 


<5sys 


bin 


range 


Born 


dressed 


bare 


/ 


{%} 


PS] 


Born 


dressed 


bare 


/ 


[%] 


[%} 


0.000 


- 0.004 


9.77 


9.69 


9.70 


1 


0.46 





35 


9.77 


9.67 


9.67 


1 


0.39 


0.28 


0.004 


- 0.008 


9.68 


9.59 


9.59 


1 


0.47 





26 


9.76 


9.66 


9.66 


1 


0.39 


0.18 


0.008 


-0.012 


9.42 


9.36 


9.38 


1 


0.47 





28 


9.42 


9.34 


9.35 


1 


0.40 


0.24 


0.012 


- 0.016 


9.14 


9.06 


9.07 


1 


0.48 





35 


9.26 


9.17 


9.18 


1 


0.40 


0.24 


0.016 


- 0.020 


8.82 


8.76 


8.77 


1 


0.49 





24 


8.83 


8.76 


8.77 


1 


0.41 


0.19 


0.020 


- 0.024 


8.48 


8.43 


8.43 


1 


0.50 





25 


8.51 


8.44 


8.45 


1 


0.42 


0.27 


0.024 


- 0.029 


7.97 


7.93 


7.94 


1 


0.46 





26 


8.05 


8.00 


8.01 


1 


0.39 


0.24 


0.029 


- 0.034 


7.57 


7.52 


7.53 


1 


0.47 





22 


7.57 


7.52 


7.53 


1 


0.40 


0.19 


0.034 


- 0.039 


7.02 


7.00 


7.01 


1 


0.49 





29 


7.11 


7.09 


7.09 


1 


0.41 


0.17 


0.039 


- 0.045 


6.55 


6.53 


6.53 


1 


0.46 





22 


6.50 


6.49 


6.49 


1 


0.39 


0.17 


0.045 


- 0.051 


5.93 


5.92 


5.92 


1 


0.48 





22 


6.00 


5.99 


5.99 


1 


0.41 


0.16 


0.051 


- 0.057 


5.52 


5.52 


5.52 


1 


0.50 





22 


5.52 


5.53 


5.53 


1 


0.42 


0.22 


0.057 


- 0.064 


5.04 


5.04 


5.04 


1 


0.48 





22 


5.00 


5.01 


5.01 


1 


0.41 


0.16 


0.064 


- 0.072 


4.55 


4.56 


4.56 


1 


0.48 





22 


4.52 


4.52 


4.52 


1 


0.41 


0.23 


0.072 


- 0.081 


4.01 


4.03 


4.03 


1 


0.48 





21 


4.04 


4.06 


4.06 


1 


0.40 


0.18 


0.081 


- 0.091 


3.58 


3.59 


3.59 


1 


0.48 





22 


3.53 


3.55 


3.55 


1 


0.41 


0.19 


0.091 


-0.102 


3.15 


3.16 


3.16 


1 


0.49 





23 


3.14 


3.16 


3.16 


1 


0.41 


0.21 


0.102 


- 0.114 


2.73 


2.74 


2.74 


1 


0.50 





26 


2.73 


2.75 


2.74 


1 


0.43 


0.23 


0.114 


-0.128 


2.34 


2.35 


2.35 


1 


0.50 





29 


2.35 


2.37 


2.37 


1 


0.42 


0.24 


0.128 


- 0.145 


2.00 


2.01 


2.01 


1 


0.49 





24 


2.00 


2.01 


2.01 


1 


0.42 


0.21 


0.145 


-0.165 


1.687 


1.697 


1.698 


1 


0.49 





28 


1.669 


1.680 


1.679 


1 


0.42 


0.28 


0.165 


- 0.189 


1.355 


1.364 


1.363 


1 


0.50 





25 


1.355 


1.365 


1.365 


1 


0.43 


0.19 


0.189 


-0.219 


1.079 


1.087 


1.087 


1 


0.50 





23 


1.086 


1.093 


1.093 


1 


0.43 


0.20 


0.219 


- 0.258 


8.27 


8.34 


8.32 


io- 1 


0.50 





24 


8.22 


8.28 


8.27 


io- 1 


0.43 


0.19 


0.258 


-0.312 


5.97 


6.00 


5.99 


io- 1 


0.50 





25 


5.94 


5.97 


5.97 


io- 1 


0.43 


0.17 


0.312 


- 0.391 


3.97 


3.99 


3.99 


io- 1 


0.51 





22 


3.94 


3.96 


3.96 


io- 1 


0.44 


0.17 


0.391 


- 0.524 


2.28 


2.28 


2.28 


io- 1 


0.52 





24 


2.29 


2.29 


2.29 


io- 1 


0.45 


0.19 


0.524 


- 0.695 


1.176 


1.179 


1.177 


io- 1 


0.64 





29 


1.164 


1.168 


1.166 


IO" 1 


0.55 


0.23 


0.695 


-0.918 


5.79 


5.80 


5.79 


IO" 2 


0.79 





37 


5.77 


5.79 


5.78 


io- 2 


0.69 


0.29 


0.918 


- 1.153 


2.94 


2.95 


2.95 


IO -2 


1.07 





47 


2.91 


2.91 


2.91 


IO -2 


0.95 


0.37 


1.153 


- 1.496 


1.54 


1.55 


1.54 


ur 2 


1.22 





52 


1.50 


1.51 


1.51 


IO -2 


1.09 


0.41 


1.496 


- 1.947 


7.25 


7.26 


7.25 


KT 3 


1.55 





66 


7.05 


7.06 


7.06 


IO" 3 


1.39 


0.48 


1.947 


- 2.522 


3.52 


3.51 


3.50 


IO" 3 


1.97 





78 


3.56 


3.56 


3.55 


IO -3 


1.74 


0.62 


2.522 


- 3.277 


1.73 


1.73 


1.72 


IO -3 


2.46 





96 


1.79 


1.80 


1.80 


10~ 3 


2.13 


0.71 



and improve the accuracy of the multi-jet background nor- 
malisation. 

The total fraction of background events is (0.61 ±0.31)% 
in the e + e~ channel and (0.56 ±0.28)% in the chan- 
nel. The multi-jet background represents ~ 50% of the to- 
tal background in both channels and dominates at low (f>* 
values. An irreducible background may also arise from the 
production of a lepton pair via photon-photon interactions, 
77 — > £ + £~. This contribution was evaluated at leading 
order using FEWZ 3.1 [24, 60] and the MRST2004qed [61] 
PDF, currently the only available PDF set containing a 
description of the QED part of the proton. According to 
the LO cross section calculated in the fiducial lepton accep- 
tance, the fraction of photon-induced events is expected to 
be below 0.1%, with an uncertainty of 50%. This contri- 
bution is six times lower than the sum of other background 
contributions and is therefore neglected. 

6. Cross-section measurement and systematic un- 
certainties 

The differential cross section is evaluated in bins of 0*, 
or of (<fi*yz), from the number of observed data events 



in each bin after subtraction of the estimated number of 
background events. 

A bin-by-bin correction is used to correct the observed 
data for detector acceptances and inefficiencies, as well as 
for QED FSR. The correction factors are determined using 
signal MC events. For the chosen bin widths the purity, 
dchned as the fraction of simulated events reconstructed 
in a (f>* bin which have generator-level (f>* in the same bin, 
is always more than 83% and reaches 98% in the highest 
4>* bins. In each bin, the data are normalised to the cross 
section integrated over the fiducial acceptance region. 

An analysis of systematic uncertainties was performed, 
in which the sensitivity of the measurements to variations 
in the efficiencies and energy scales of the detector com- 
ponents and to the details of the correction procedure is 
tested. The systematic uncertainties in the measured cross 
section are determined by repeating the analysis after ap- 
plying appropriate variations for each source of systematic 
uncertainty to the simulated samples. The systematic un- 
certainties which are correlated between (j)* bins are listed 
below. 

• Uncertainties in the estimation of the number of back- 



4 



ground events from multi-jet, W — > Iv and — > 
t + t~ decays, tt and diboson processes yield values 
of up to 0.3% in the e + e~ and ^t + yU _ channels, when 
propagated to the normalised differential cross sec- 
tion. 

• Possible mis-modelling of the angular resolution of 
tracking detectors leads to uncertainties of up to 0.3% 
(0.2%) on the normalised differential cross section in 
the e + e~ (fi + jX~ ) channel. 

• The dependence of the bin-by-bin correction factors 
on the shape of the assumed </>* distribution was 
tested by re-weighting simulated events to the mea- 
sured </>* cross section. An iterative Bayesian unfold- 
ing technique [62] was employed as an alternative ap- 
proach to assess systematic uncertainties. The un- 
certainty in the correction procedure is found to be 
smaller than 0.1% in both channels and for the full 
(f)* range. 

• As the definition of the 0* variable is based on the 
lepton angles, the normalised differential cross sec- 
tion depends only weakly on uncertainties in the lep- 
ton energy/momentum scale and resolution. When 
propagated to the normalised differential cross sec- 
tion, these uncertainties amount to less than 0.1% and 
0.03% in the e + e~ and channels, respectively. 

• Uncertainties arising from the mis-modelling of lepton 
identification efficiencies and trigger efficiencies in the 
simulation amount respectively to 0.05% (0.03%) and 
0.04% (0.02%) in the e+e" channel. 

• Pilc-up has only a weak influence on this measurement 
and results in an uncertainty of at most 0.05% on the 
normalised differential cross section. 

A second class of systematic uncertainties, listed below, 
are considered uncorrelated across </>* bins. 

• Uncertainties on the bin-by-bin correction factors 
arising from the MC sample statistics are 0.2% 
(0.13%) at low (j>* in the e + e~ (fi + [j,~) channel, in- 
creasing to 0.9% (0.6%) in the highest cf>* bins. 

• Possible local biases in angular measurements (</>, rj) 
by tracking detectors yield an estimated constant un- 
certainty of 0.1% on the normalised differential cross 
section. The local effect of these biases allows bin-to- 
bin correlations to be neglected. The impact of this 
assumption on the combination of electron and muon 
channel results is small. 

• A conservative systematic uncertainty of 0.3% due 
to c/>*-dcpendcnt modelling of QED FSR is assigned 
by comparing predictions from Photos [49] and 
from the Sherpa implementation of the YFS algo- 
rithm [50, 51]. This comparison provides the size 
of the uncertainty but however does not allow the 



shape of the </>* dependence to be estimated. The 
uncertainty is assumed to hold for cross sections at 
Born, dressed and bare levels and for both electron 
and muon channel measurements. It therefore does 
not affect the combination of them. 

The total systematic uncertainty on each data point is 
formed by adding the individual contributions in quadra- 
ture. 
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Figure 1: The measured normalised differential cross section l/<T fld • 
dcr fid /d0* as a function of <f>* for Z/'y* — > e+e~ (closed dots) and 
Z/~/* — ¥ /x+ fi~ (open dots) channels. The measurements are com- 
pared to ResBos predictions represented by a line. The ratio of 
measured cross sections to ResBos predictions is presented in the 
bottom panel. The measurements are displaced horizontally for bet- 
ter visibility. The inner and outer error bars on the data points 
represent the statistical and total uncertainties, respectively. The 
uncertainty due to QED FSR is included in the total uncertainties. 



7. Results and discussion 



The normalised differential cross sections measured for 
Z/^f* — > e + e~ and Z/"f* — > /i + /-t _ production in the fidu- 
cial acceptance are presented in Table 1. The measure- 
ments are reported with respect to the Born, dressed and 
bare reference points at particle level regarding QED FSR. 
The QED FSR corrections for the three levels are calcu- 
lated using Photos. The measured cross sections defined 
at the Z/"f* Born level are shown in Fig. 1 for the e + e _ 
and fi + jX~ channels and are compared to predictions from 
ResBos. 
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The normalised differential cross sections measured in 
the fiducial acceptance for the two channels are combined 
using a % 2 minimisation method which takes into account 
the point-to-point correlated and uncorrelated systematic 
uncertainties [63-65] and correlations between electron 
and muon channels. The uncertainties due to the unfold- 



Table 2: The combined normalised differential cross section l/a^ d ■ 
do- fld /d</>* in bins of <f>* at Born level. The statistical (<5 s tat) and total 
systematic (<5 sys ) uncertainties are given in percent. The normalised 
differential cross section extrapolated to the full lepton acceptance 
l/o- tot -dir t ° t /d(/>* is obtained at Born level by multiplication with the 
inverse acceptance correction factor A^ 1 ■ The uncertainty <5(AJ ) 
on this acceptance correction factor is also given in percent. The 
overall point-to-point uncorrelated additional uncertainty in QED 
FSR of 0.3% is not included. 
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interfaced to Pythia8. Uncertainties in A c resulting from 
PDF uncertainties are below 1%. 
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Figure 2: The ratio of the combined normalised differential cross 
section l/<r tid ■ do- fld /d0* to ResBos predictions as a function of <f>*. 
The inner and outer error bars on the data points represent the sta- 
tistical and total uncertainties, respectively. The uncertainty due to 
QED FSR is included in the total uncertainties. The measurements 
are also compared to predictions, which are represented by a dashed 
line, from Ref. [22] and from Fewz in the top and bottom panels, 
respectively. Uncertainties associated with these two calculations 
are represented by shaded bands. The prediction from Fewz is only 
presented for <f>* > 0.1. 



ing procedure, the pile-up, and QED FSR arc considered 
to be completely correlated between the e + e~ and /i + /-t~ 
channels. The minimisation yields a total x 2 per degree 
of freedom (ridof) of x 2 / n doi = 33.2/34, indicating a good 
consistency between the electron and muon data. Mea- 
sured values of the combined normalised differential cross 
section l/<r fid • dcr fid /d0* within the fiducial lepton accep- 
tance are presented in Table 2. The acceptance correction 
factors A c needed to extrapolate the measurement to the 
full lepton acceptance are determined using Powheg and 
the CT10 PDF set. The uncertainty in A c is estimated 
from the differences among predictions obtained with Res- 
Bos, Mc@nlo, Sherpa, Alpgen, Herwig and Powheg 



The ratio of the combined normalised differential cross 
section to the ResBos prediction is shown as a function 
of (f>* in Fig. 2. The measurement is also compared to a 
QCD calculation by A. Banfi et al. [22] and to another 
obtained with Fewz 2.1. The ratios of these two calcu- 
lations to ResBos predictions are also shown in Fig. 2. 
The CTEQ6m [66] PDF set is used in the calculations in 
Ref. [22] . The theoretical uncertainties on this calculation 
are evaluated by varying the resummation, renormalisa- 
tion and factorisation scales /^q, [1r and \if between mz/2 
and 2mz, with the constraints 0.5 < fAi/fJ-j < 2, where 
i,j € {F,Q,R} 7 and /if/^q > 1- Uncertainties coming 
from the PDFs are also considered [22]. For Fewz, the 
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Table 3: The combined normalised differential cross section l/cr fld ■ d<r fld /d0* in bins of </>* and in three \yz\ ranges. The statistical (<5 s tat) 
and total systematic (S ByB ) uncertainties are given in percent. The overall point-to-point uncorrelated additional uncertainty in QED FSR of 
0.3% is not included. 



bin 



range 



\yz 



< 0.8 



l/cr hd ■ dcr hd /dcf>* 



<5stat 
[%] 



"svs 
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0.8 < \y z \ 
l/o- fid • der fld /d0* 



< 1.6 

<5stat 
[%] 



[%] 



|yz| > 1.6 

l/ CT fld -do- fld /d<#>; 



<5stat 
[%] 



"svs 

[%] 



0.000 

0.004 

0.008 
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0.016 
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0.034 
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0.391 

0.524 
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1.496 

1.947 - 

2.522 
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CT10 PDF set is used. Uncertainties are evaluated by 
varying and by factors of two around the nominal 
scale mz with the constraint 0.5 < \irI [If < 2, by varying 
a s within a range corresponding to 90% confidence-level 
(CL) limits [67], and by using the PDF error eigenvector 
sets. 

The difference between the ResBos prediction and data 
is ~ 2% for (f>* < 0.1, increasing to 5% for higher </>* values. 
This difference is smaller than the uncertainty in ResBos 
predictions due to the propagation of PDF eigenvectors 
sets, which amounts to 4% for </>* < 0.1 and 6% above. The 
description of data provided by calculations from A. Banfi 
et al. [22] is less good than ResBos but observed differ- 
ences remain within the theoretical uncertainties of the 
calculation. The prediction obtained with Fewz under- 
shoots the data by ~ 10%, as already observed for the p^ 
spectrum in Ref. [15]. At low (f>* values, corresponding 
mainly to low p^, fixed-order perturbative QCD calcula- 
tions are not expected to give an adequate description of 
the cross section. The prediction from Fewz is therefore 
only presented for <f>* > 0.1. It is normalised using the 
total cross section predicted by Fewz, which accurately 
describes experimental measurements [58]. 



The cross section is also measured double differentially 
in bins of (f>* for three independent bins of \yz\ for both the 
e + e~ and fi + fi~ channels. The double differential cross- 
section measurements in the two channels are combined 
using the same x 2 minimisation procedure as used for the 
single differential cross section. The minimisation yields 
a total x 2 / n dof = 118/102. Measured values of the com- 
bined normalised differential cross section l/<7 fid -dcr fid /d</>* 
within the fiducial lepton acceptance in all (/)* and \yz \ bins 
are presented in Table 3. 

The ratio of the combined normalised differential cross 
section to the ResBos prediction is shown as a function 
of </>* for the three \yz\ ranges in Fig. 3. The measure- 
ment is also compared to predictions obtained using dif- 
ferent MC event generators. The PDF set CT10 is em- 
ployed in all calculations, except for Alpgen where the 
CTEQ6L1 PDF set is used. The parton-shower parame- 
ters of each MC generator are set to their default values, 
except for Pythia6 where a specific ATLAS re-tuning was 
used [48]. The generators Alpgen, interfaced to Herwig, 
and Sherpa provide a good description of the spectrum 
for 4>* > 0.1. In particular, Sherpa describes the data 
better than ResBos over all \y z \ bins for 4>* > 0.1. How- 
ever, for 0* < 0.1 the deviations of Sherpa or Alpgen 
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Figure 3: The ratio of the combined normalised differential cross section l/cr fld ■ d<r fid /d</)* to the ResBos predictions as a function of (j>* in 
three ranges of \yz\- The inner and outer error bars on the data points represent the statistical and total uncertainties, respectively. The 
uncertainty due to QED FSR is included in the total uncertainties. The measurements are also compared to predictions from different MC 
event generators. 



from the data are ~ 5%, somewhat larger than those of 
ResBos. The Powheg generator interfaced to Pythia8 
is also able to describe the data to within 5% over the 
whole 4>* range. 

The effect of changing the PS tunings and algorithms in- 
terfaced to Powheg was investigated by using Pythia6 
and HERWIG interfaced to the same Powheg NLO cal- 
culation. These two variations give a worse description of 
data than Pythia8, and deviations from data of ~ 10% 
are observed. The Mc@nlo generator interfaced to Her- 
WIG does not properly describe the data for <jj* > 0.1, and 
deviations from data of the order of 4-7% are observed for 
4>* < 0.1 depending on the \yz\ bin. The level of agree- 
ment between MC generators and data is very similar for 
comparisons at the dressed level. 



8. Conclusion 

A measurement of the </>* distribution of Z/j* boson 
candidates in y/s = 7 TeV pp collisions at the LHC is pre- 
sented. The data were collected with the ATLAS detector 
and correspond to an integrated luminosity of 4.6 fb -1 . 



Normalised differential cross sections as a function of 0* 
have been measured in bins of the Z boson rapidity yz up 
to (f>* ~ 3 for electron and muon pairs with an invariant 
mass 66 GeV < m« < 116 GeV. The high number of Z/ 7 * 
boson candidates recorded permits the use of finer bins as 
compared to a similar study performed at the Tevatron. 
The typical uncertainty achieved by the combination of 
electron and muon data integrated over the whole Z ra- 
pidity range is below 0.5% for 0* < 0.5 increasing to 0.8% 
at larger <p* values. 

The cross-section measurements have been compared 
to resummed QCD predictions combined with fixed-order 
perturbative QCD calculations. Calculations using Res- 
Bos provide the best descriptions of the data. However, 
they are unable to reproduce the detailed shape of the 
measured cross section to better than 4%. 

The cross-section measurements have also been com- 
pared to predictions from different Monte Carlo genera- 
tors interfaced to a parton shower algorithm. The best 
descriptions of the measured <p* spectrum are provided by 
Sherpa and Powheg +Pythia8 Monte Carlo event gen- 
erators. For <f>* values above 0.1, predictions from Sherpa 



are able to reproduce the data to within ~ 2%. The low cj>* 
part of the spectrum is, however, described less accurately 
than by ResBos. Double differential measurements as a 
function of 0* and yz provide valuable information for the 
tuning of MC generators. None of the tested predictions is 
able to reproduce the detailed shape of the measured cross 
section within the experimental precision reached, which 
is typically lower by one order of magnitude than present 
theoretical uncertainties. 
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